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Charge transfer in collisions of Na^ cluster ions with Cs 
atoms is investigated theoretically in the microscopic frame- 
work of non-adiabatic quantum molecular dynamics. The 
competing reaction channels and related processes affect- 
ing the charge transfer (electronic excitations, fragmentation, 
temperature) are described. Absolute charge transfer cross 
sections for Na^(2.7 keV) -I- Cs — > Na^ -I- Cs"*" have been 
calculated in the size range 4 < n < 11 reproducing the size 
dependence of the experimental cross sections. The energy 
dependence of the cross section is predicted for n = 4, 7, 9. 
An exotic charge transfer channel producing Cs~ is found to 
have a finite probability. 

34.70.-He, 36.40.Qv, 31.15.Qg 



Charge transfer represents one of the fundamental 
atomic interactions. In the last decades, large progress 
has been made in the understanding of charge transfer 
mainly for two borderline cases: the elementary ion~atom 
reaction ^ and the complex ion-surface interaction 
In the recent past, however, there has been a flurry of 
activity to close the gap between the two fields by in- 
vestigating the intermediate case of ion-cluster collisions 
|p|-|l8|. In general, the basic aspect of cluster collisions 
consists in the simultaneous occurrence and mutual cou- 
pling of electronic transitions (charge transfer, excitation, 
ionization) and nuclear-core excitations (vibration, rota- 
tion, fragmentation) in a system, where the number of 
electronic and nuclear degrees of freedom is large hut fi- 
nite. In particular, collision induced dissociation (CID) 
will compete with charge transfer (CT) in such collisions. 

Already, in one of the first experiments with mass- 
selected cluster beams, CT and CID were investigated 
for Na+ + Cs |,| and K+ -f Cs [| collisions. An im- 
portant result of these novel experiments was the direct 
determination of the CT cross section for the neutraliza- 
tion of the parent cluster ion, i.e. 



Na^ 



Cs 



Na„ + Cs^ 



(1) 



by measuring the neutral products associated with the 
process (^) but discriminating neutral products originat- 
ing from CID. Further experimental studies of CT in col- 
lisions of cluster cations with atoms ^-p^ , molecules |^] 



and clusters pl|,p[ have been performed providing a con- 
siderable amount of data. Multiple CT processes accom- 
panied by fragmentation have been observed in collisions 
of highly charged ions with fuUerenes and sodium 

clusters Q. 

This enormous experimental progress was accompa- 
nied by the development of a number of theoretical 
descriptions of CT reactions in cluster collisions. In 
Ref. Q, a two-state model of (near) resonant CT was 
presented. Classical barrier models ||l^,|l^ have been ap- 
plied to distant collisions of Ceo with highly charged ions. 
Semi-microscopic descriptions of CT, which are based on 
the jellium approximation [|l7| or on phenomenological 



single-particle potentials [IS), have also been suggested. 
However, the combined description of CT and fragmenta- 
tion or, more generally, a simultaneous treatment of elec- 
tronic and nuclear degrees of freedom in non-adiabatic 
cluster collisions is still an unsolved problem. 

Recently, a universal microscopic approach called non- 
adiabatic quantum molecular dynamics (NA-QMD) has 
been developed, which describes classical atomic motion 
simultaneously and self-consistently with electronic tran- 
sitions in atomic many-body systems ]T9| . In this pa- 
per, we present a fully microscopic analysis of CT and 
fragmentation in cluster collisions applying the NA-QMD 
theory to Na+ -f Cs. 

The NA-QMD approach has been derived [|l9| on the 
basis of the time-dependent density functional theory 
pO| , whereby the Kohn-Sham formalism within the time- 
dependent local density approximation is used. Resulting 
from an LCAO ansatz for the Kohn-Sham orbitals a set 
of coupled differential equations for the time-dependent 
coefficients is obtained to determine the time evolution 
of the electronic density as the consequence of the classi- 
cal atomic motion. Simultaneously, Newton's equations 
of motion with explicitly time-dependent forces have to 
be solved reflecting the possible energy transfer between 
the classical system of ionic cores and the quantum- 
mechanical system of valence electrons. Details of the 
procedure to calculate expectation values of observables 
or probabilities for specific electronic transitions (defined 
by the time-dependent electronic density) as well as of the 
different fragmentation channels (obtained from classical 
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trajectories) will be given elsewhere 

To elucidate the different physical processes determin- 
ing the absolute cross section of the CT reaction (|^), 
we have performed a detailed study of the system NaJ 
+ Cs at a lab collision energy of E^h = 2.7 keV. The 
cluster projectile is prepared initially in its electronic 
and geometric ground state (rhombic, D2h)- In Fig- 0, 
the time evolution of the calculated mean charge lo- 
cated at the Cs atom {qcs){t) (obtained from a popu- 
lation analysis) as well as of the kinetic-energy differ- 
ence AEi.init) = -Ec.m. -£'kin(i), whcrc Ec.m. and Ei^init) 
are the collision energy and the total kinetic energy, re- 
spectively, both referred to the center-of-mass frame, are 
presented for two collision events with the same impact 
parameter (b = 9a.u.) but different initial orientations of 
the cluster with respect to the beam axis (see inserts (a) 
and (b) in Fig. |l|). The quantity {qcs){t) characterizes the 
quantum dynamics of CT, whereas Ai?kin(0 describes 
the energy flow in the classical (ionic) degrees of free- 
dom. A strong dependence of both quantities on the ini- 
tial orientation of the cluster is observed. In the example 
(a), nearly compensating charge fluctuations during the 
interaction lead finally to (qcs) ~ 0, and the event corre- 
sponds to a nearly elastic scattering AE'kin « 0. In con- 
trast, the large CT in the case (b) (qcs) ~ 0.75 is accom- 
panied by a considerable energy loss of AiJkin ~ 0.34 eV. 
An endothermic character is typical for most collision 
events with appreciable CT, even at larger impact pa- 
rameters, where no vibrational excitation of the cluster 
occurs. Comparing the ionization potentials of the Cs 
atom (3.89 eV) and the Na4 cluster (4.24 eV ^]) and 
assuming a CT that results in the electronic ground state 
of Na4, one should expect an exothermic reaction with 
AE'kin ~ —0.35 eV. That this is not the case indicates 
that the CT process produces an electronically excited 
Na4 cluster. 

The calculation of CT cross sections requires a de- 
tailed analysis of the final electronic and atomic states 
and, in particular, the careful consideration of CT prob- 
abilities for different reaction channels. To start with a 
transparent classification, we define integral CT proba- 
bilities, which are the probabilities P(Cs'^) to find the Cs 
atom in the charge state q in the exit channel, where q 
is an integer and J2q P{Gs'^) = 1. These probabilities 
describe the primary CT (q ^ 0) and scattering {q = 0) 
processes without regard to the further evolution of the 
cluster, i.e. to the possible fragmentation. The related 
CT cross sections can be directly measured by detect- 
ing the formed Cs ions. The probabilities P(Cs'^), cal- 
culated as an average of the results obtained with about 
300 different initial orientations of the cluster per im- 
pact parameter b, are shown as a function of b in the 
upper panel of Fig. ^. The CT leading to Cs"^ ions and 
the scattering without CT have nearly equal probabilities 
P(Cs~'') and P(Cs), respectively, for impact parameters 
6 < 8 a.u.. As the graph of P(Cs^) indicates, CT takes 



place with remarkable probability up to impact param- 
eters of 6 « 15 a.u., which is more than twice the long 
half-axis {R = 5.7a.u.) of the Na4 rhombus. Surprisingly, 
the calculations also yield a finite probability P(Cs~) for 
an electron transfer to the Cs atom, which represents an 
interesting prediction for future experimental studies. 



Na/ (2.7 keV) + Cs (b=9a.u.) 
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FIG. 1. Calculated kinetic-energy difference ABkin (upper 
panel) and mean value of the charge located at the Cs atom 
(gcs) (lower panel) as functions of time t for the two initial 
collision geometries illustrated in the inserts (a) and (b). The 
dashed and solid lines correspond to the cases (a) and (b), 
respectively. 

The 6-weighted integration of P(Cs^) leads to a to- 
tal CT cross section of cr(Cs+) = (38.2 ± 2.1) For 
comparison, the "geometrical" cross section of the clus- 
ter is about do ~ 71" P^ « 29 A and the measured value 
is cr°''P(Na4) = (17 ± 3) A^ |]. It is important to re- 
alize, however, that, instead of Cs+, the signal of the 
neutralized Na4 was detected in the measurements 
The "survival" probability of the Na4 cluster after the 
primary CT is given by P(Cs^) • [1 — Ppr], where Ppr is 
the total fragmentation probability. In cluster collisions, 
fragmentation can be induced through three mechanisms 
characterized by different time scales [ ^l|j2^ , p^ . These 
are large momentum transfer between atoms of the pro- 
jectile and of the target, electronic excitation followed 
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by energy transfer via electron- vibrational coupling, and 
statistical fragmentation. Whereas the first two mech- 
anisms are precisely described by the NA-QMD theory 
and automatically accounted for in the actual calcula- 
tions (this has been shown in comparisons with experi- 
mental data on CID ||2^), the total probabilities of the 
different fragmentation channels can be determined from 
the calculated internal (vibrational and electronic) exci- 
tation energy of the cluster using statistical arguments 
pl[ . The calculated total fragmentation probability Ppr 
and the probability P(Na4) to have a neutral Na4 cluster 
in the exit channel are shown as a function of b in the 
lower panel of Fig. 
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FIG. 2. Calculated integral CT probabilities P(Cs+), 
P(Cs), P(Cs~) [upper panel], and fragmentation Ppr and CT 
probabilities P(Na4) [lower panel] as functions of the impact 
parameter b for NaJ chosen initially at T = K in its low- 
est-energy structure (see text). The length of the error bars 
is given by 2 so, where so is the standard deviation of the 
orientation average. 

Fragmentation occurs with high and then decreasing 
probability up to 6 « 11 a.u., and the computed to- 
tal fragmentation cross section of ap^ = 49.4 A ex- 
ceeds the geometrical cross section considerably. The 
CT probability P(Na4) peaks around b ^ 9 a.u., which is 
larger than the cluster "radius" R. The CT cross section 
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CT(Na4) = 20.2 A is still slightly larger than the experi- 
mental value. 

The calculations presented so far have been carried out 
with zero initial temperature of the cluster, whereas un- 
der the experimental conditions |^,^ "liquid" cluster ions 
are used resulting from the laser ionization to produce 
the cluster beam |25 . To take into account the temper- 



ature effect, we have repeated the procedure described 
above with excited clusters in the initial state, where the 
total excitation energy was chosen to be slightly below 
the dissociation energies of the neutral as well as of the 
cationic cluster. The initial configurations of the excited 
clusters for the simulation of the collisions were produced 
in equilibration runs of 300 ps in length. The calculated 
CT cross section cr(Na4) — 16. SA obtained with "liquid" 
cluster ions is in perfect agreement with the experimental 
result. 

The theoretical and experimental j|] CT cross sections 
are compared as a function of the cluster size n in Fig. 
The theoretical results were obtained with "liquid" clus- 
ter ions in the initial state. Except for the case of n = 5, 
the agreement of the computed and measured data can 
be qualified as perfect, since these data are absolute cross 
sections. The statistical uncertainties of the calculated 
cross sections reflect the strong dependence of the out- 
come of the collision process on the initial configuration. 
These uncertainties could be reduced further only at a 
very high computational cost The overestimated 

cross section in the case n — 5 results from contributions 
of a particular isomer, for which accidentally a (quasi-) 
resonant CT with Cs occurs 

Na,"(2.7 keV) + Cs ^ Na, + Cs' 




n 

FIG. 3. Calculated and measured [^ CT cross sections 
CTCT(Na„). The theoretical error bars result from the sta- 
tistical uncertainty 2 sd (cf. Fig. ^ . 

The distinct maximum of the experimental CT cross 
section at n = 7 is reproduced by the calculations. This 
maximum, however, may be peculiar to the particular 
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collision energy. In the experimental investigation of 
K+ + Cs collisions [Q a very strong dependence of the 
CT cross section on the collision energy has been found 
(cf. Fig. 4 in Ref. [Q). In order to examine this aspect 
and to stimulate further experimental investigations, we 
have calculated the CT cross section over a wide range 
of the (center-of-mass) collision energy (0.1 ... 30 keV) 
for n = 4, 7, 9. The results are shown in Fig. ^. The 
absolute CT cross section for n = 7 exceeds those for 
n = 4 and n — 9 over the entire energy range considered. 
Consequently, the large CT cross section in Na^^ + Cs 
collisions should be attributed to the specific electronic 
structure of Na^ providing favorable, i.e. near resonant 
conditions for CT in collisions with Cs. 



Na " + Cs ^ Na„ + Cs' 
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FIG. 4. Calculated collision energy dependence of the CT 
cross section crcT(Na„) for n = 4,7,9 (full symbols). The 
experimental data (open symbols with error bars; [^]) corre- 
spond to the same lab collision energy of iJiab ~ 2.7 keV. 

In summary, we have presented results of a fully micro- 
scopic analysis of CT and fragmentation in cluster-ion - 
atom collisions based on NA-QMD simulations. The de- 
tailed study of Na^(2.7 kcV) + Cs collisions revealed the 
role of the different physical processes associated with 
CT in cluster collisions (electronic excitations, fragmen- 
tation, temperature effects) . An exotic "inverse" CT pro- 
cess leading to Cs~ is predicted. The calculated absolute 
CT cross sections for Na+(2.7 keV) + Cs {4 < n < 11) 
are in good agreement with the experimental data. The 
energy dependence of the CT cross section is predicted 
for several cluster sizes {n — 4, 7, 9) in order to encourage 
further experimental studies. 
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